REPORT 972 


THE EFFECT OF TORSIONAL FLEXIBILITY ON THE ROLLING CHARACTERISTICS 
AT SUPERSONIC SPEEDS OF TAPERED UNSWEPT WINGS 

By Warren A. Tuck.ee and Robert L. Nelson 


SUMMARY 

An analysis is presented of the effect of torsional flexibility 
on the rolling characteristics at supersonic speeds of tapered 
unswept wings with partial-span constant-percent-chord ailer- 
ons extending inboard from the wing tip. The geometric var- 
iables considered are aspect ratio, taper ratio, aileron span, and 
aileron chord. 

The shape of the wing-torsional-stiffness curve is assumed 
and the twisting moment is considered to result solely from the 
pressure distribution caused by aileron deflection, so that the 
necessity of using a successive-approximation method is avoided. 

Because of the complexity of the equations resulting from the 
analysis, numerical calculations from the equations are pre- 
sented in a series of figures. A computational form is pro- 
vided to be used in conjunction with these figures so that 
calculations can be made without reference to the analysis. 

INTRODUCTION 

In reference 1 an analysis was presented of the effect- of 
torsional flexibility on the rolling characteristics at super- 
sonic speeds of rectangular wing-aileron combinations. 
The present report deals with the more general case of 
linearly tapered wings having straight midchord lines and 
incorporating partial-span constant-percent-chord ailerons 
extending inboard from the wing tip (fig. 1) and may thus 
he considered a companion piece to reference 2, which gives 
an analysis of the rolling characteristics of this class of 
wings at subsonic speeds. 

The necessity of using successive approximations is avoided 
by assuming, as in reference 1, that the principal twisting 
moment is contributed by the loads arising from aileron 
deflection. The various aerodynamic forces and moments 
are calculated by methods based mainly on the linearized 
equations for three-dimensional supersonic flow. 

The geometric variables considered are aspect ratio, taper 
ratio, aileron span, and aileron chord. The torsional stiff- 
ness of the wing is assumed to vary inversely with the cube 
of the distance outboard from the wing center line. 

Because of the length of the equations for the various 
quantities appearing in the analysis, numerical calculations 
have been made for ranges of the variables considered. 
The results can be used, together with a computational 
form that is provided, to make calculations without refer- 
ence to the details of the analysis. 
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SYMBOLS 


aspect ratio 



b* 2b 

S Cr(l-rX) 


) 


wing span 

total span of two ailerons 

wing chord at spanwise station y 

aileron chord at spanwise station y 

root chord 

tip chord 

taper ratio (c ( /c r ) 


slope of hinge line (see fig. 1) ^ 


slope of leading edge (see fig. 1) 



wing area 



semispan (5/2) 
coordinate axes (see fig. 1) 


M . 

My 


free-stream klach number 
free-stream Mach number component normal 
to hinge line 
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e r 

pb\2V 

p 


w 

V 

L 

Oi 

ft, 


ft 


V 


ratio of specific heat at constant pressure to 
■ specific heat at constant volume (1.40 
for air) 

free-stream static pressure 


free-stream dynamic pressure 



lifting pressure 

lifting-pressure coefficient (P/q) 
additional two-dimensional lifting-pressure 
coefficient due to aileron deflection 



deflection of one aileron, measured in plane 
perpendicular to wing and parallel to free 
stream, positive when right aileron is 
deflected up 

angle of twist at any spanwise station, 
positive when leading edge of right wing 
moves up 

angle of twist at reference station (midspan 

of aileron) 

wing-tip helix angle 

rolling angular velocity, positive when right 
wing moves down 

constant of proportionality for parabolic 
twist distribution 
vertical disturbance velocity 
free-stream velocity 

rolling moment, positive when tending to 
depress right wing 
rolling-moment coefficient ( L/qSb ) 
rolling-moment-effectiveness coefficient, for 
two ailerons (jbCifbS) 


damping-in-roll coefficient 




T'=T 


2/3 

Op a qc t 3 


rolling-moment-loss coefficient (d(7j/50 r ) 
aerodynamic twisting moment at any span- 

wise station ( r, md *) 


m pitching moment about midchord at any 

spanwise station 

6J torsional-stiffness parameter 



me torsional-stiffness parameter (T/d) 

me T torsional-stiffness parameter at reference 

station (midspan of aileron) 

k constant of proportionality for stiffness 

distribution 

F trailing-edge-angle reduction factor 

4 > included trailing-edge angle measured parallel 

to free stream 

fa included trailing-edge angle measured normal 

to hinge line 


Subscript: 

R rigid 

0 sea level 

1, 2, 3, 4 indicate various regions of the wing 

a, b, o indicate various stations 

All angles are in radians, unless otherwise specified. 

ANALYSIS 

GENERAL 

In a. steady roll, the resultant rolling moment acting on a 
wing of any plan form is zero. This statement can be written 
in equation form as : 

f3Ci t 8+j3C tp !£+j8<7, #r * r -0 (1) 

Equation (1) can be solved either for the rate of roll per 
unit aileron deflection or for the ratio of the rate of roll of a 
flexible wing to the rate of roll of a perfectly rigid wing 
(0 r =O) for a given aileron deflection. The following equa- 
tions are obtained: 


pb/2V P C h+P C h T l 
5 =pUT r 

pb/2V -P G % 9 r 

(p6/ 2V) B f3Ci s 8 


( 2 ) 

(3) 


The problem now becomes one of finding expressions for the 
quantities pC is , /3 C lp , / 3C le , and 8 r /8 for the class of wing- 

aileron combination considered; the rest of the analysis is 
concerned with this problem. 

A knowledge of the twist along the span is needed for the 
determination of pC, a and 0 r /5; this requirement in turn 

necessitates finding the aerodynamic twisting moment at any 
spanwise- station. . In order to evaluate the aerodynamic 
twisting moment the following assumptions are made : 

(1) The elastic axis and the center of pressure at any 
spanwise' station of the loads resulting from rolling and from 
twist lie very near the midchord line. This assumption 
allows the twisting moments resulting from rolling and from 
twist to be neglected. 

(2) The aerodynamic twisting moment results solely from 
the loads caused by aileron deflection. This assumption is a 
direct result of the first one. 

(3) The loads caused by aileron deflection, wing twist, and 
rolling are mutually independent. This assumption corre- 
sponds to the familiar superposition principle of linearized 
theory. 

(4) The aileron does not twist with respect to the wing. 

An important result of making these assumptions is the 

elimination of any necessity for using a method of successive 
approximations. 

The main part of the analysis follows: First, the pertinent 
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pressure distributions are discussed; nest, the method of 
finding the twist along the span is outlined; and last, the 
results for the various quantities are given. Since the 
mathematics involved in each step is lengthy, no attempt is 
made to give all the details. Those readers whose primary 
interest is in the application of the results are advised to 
turn to the section entitled “Computational Procedure and 
Discussion” and to table I. 


PRESSURE DISTRIBUTIONS 

Three pressure distributions are necessary in order to de- 
termine jS Oi s , and pC i} . These pressure distributions 

are, respectively, that due to aileron deflection, that due to 
a wing twist increasing linearly with distance along the span 
(corresponding to a rolling wing), and that due to a wing 
twist increasing with the square of the distance along the 
span (the reason for using a parabolic twist is discussed in 
the section entitled “Twist ratio”). Each pressure distri- 
bution can be found by use of the methods of references 3 
and 4 if appropriate source-sink distributions are used to 
represent in one case the deflected aileron and in the other 
cases the twisted wing. The details of the derivations are 
not given; an illustration of the use of the method for finding 
the pressure distribution over a twisted wing can be found 
in the appendix of reference 1 . 

The equations for the pressure distribution due to aileron 
deflection are 



for the regions indicated in sketch 1. 


(4) 



twist increasing linearly with distance along the span 
(w=py) are 
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where 


Py 

x 


v»= 


( s—y ) 


x—k«s 

and the regions are indicated in sketch 2. 


(5) 



The equations for the pressure distribution due to a wing 
twist increasing with the square of the distance along the 
span ( w—ay 2 ) are 


[( 2 + ? + 5 1 ) cos “‘ a/S 
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where vi and v 2 and the regions (sketch 2) are the same as for the wing with twist increasing linearly with distance along 
the span. 


TWIST 


Aerodynamic twisting moment. — The aerodynamic twisting moment T at any spanwise station y is simply the total 
moment acting outboard of the station. Thus, 

/>• 

(7) 


T = J m iy 


where m, the pitching moment about the midchord line at any spanwise station, arises from the pressure distribution 
caused by aileron deflection and can be calculated from equations (4). 

The resulting expressions for the aerodynamic twisting moment were extremely complicated. The twisting-moment 
distribution was therefore approximated by a series of straight lines as shown in figure 2, where the exact and approxi- 
mate distributions for a particular configuration are compared for two values of j8. The equations for T' a , T' b , and T' c (see 
fig. 2) are as follows: 
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The spanwise stations for T' a , T\, and T' c are given by the following equations: 
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1 
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(9) 


1-X l-£ c 1-| b 

The straight-line approximations to the twisting-moment distribution (see fig. 2) are given by: 

T\=T' e 
T' 2 =T' e - 


T\=T\ 
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(a) 0=1.0. (b) 0=0.387. 

Figcee 2.— Comparison of esaet and approximate twisting-moment distribution tor a particular configuration. X=0.S; .4. =4; 


Twist at any spanwise station. — The twist 6 at any span- 
wise station y can be calculated from the equation 


’“If ^ 


where ddjdy can be found from the definition of the torsional- 
stiffness parameter GJ, which is 

0J -im < 12 > 

and which is related to the torsional-stiffness parameter me 
by the following equation: 

m = fo W (13) 

If the assumption is made that 

-7Tj=ky 2 (14a) 


1 k , 


This variation of torsional stiffness is known to be typical 
of many subsonic fighter airplanes (reference 2) and may be 
expected to approximate the stiffness distribution of airplanes 
or missiles designed for supersonic flight. 

If equation (14a) is substituted in equation (12) and the 
result substituted in equation (11), then the twist is expressed 
as 

d(y)=kj* y*T(y) dy (15) 

If a change is made in the variable of integration by writing 


and if T" is defined as 

2/3 

rpi rp f_L 
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then the expression for 8 becomes 

„ s C,JqsW C” an 








430 


REPORT 972 — NATIONAL ADVISORY COMMITTEE EOR AERONAUTICS 


The integration must of course be carried out in steps, since 
T' is not defined continuously over the span. 

The subscripts 1, 2, 3, and 4 may conveniently be used to 
indicate the regions designated by the same subscripts in 
figure 2. After substitution of the values of T' from equa- 
tions (10) in equation (17) the twist equations become: 
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Twist ratio. — For purposes of calculation, the twist may be 
■expressed as the ratio of the twist 6 at any station to the 
twist B r at a reference station taken at the midspan of the 

aileron (y = 1 — ^ The twist is plotted in this manner in 

figure 3 for a particular configuration for two values of d- 
The exact results are based upon the approximate twisting- 
moment distribution shown in figure 2. Because 6J9 r is not 
defined continuously over the span and because the defining 
expression is very complicated in each region, the subse- 
quent evaluation of dO s is greatly simplified if a simple 

■approximate equation for 0/0 r , continuous over the span, is 
used. This approximation, which is also plotted in figure 3, is given by the following equation: 

6 1 
6r 



Figure .3.— Comparison of exact and approximate (equation (19)) twist distribution for a 
particular configuration* \-0.5; ^-0.5; — -0.2. 
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DERIVATIVES 

Rolling-moment-effectiveness, coefficient dCy — The rolling-moment-effectiveness coefficient (3G ls is first found for a 

wing-aileron configuration of zero thickness (or zero trailing-edge angle) by suitable integration of the pressure distribu- 
tion given by*;equations (4). A factor F is then applied to account for the effect of the trailing-edge angle. 

The result obtained for the wing of zero thickness is given by the following equation: 
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The factor F that appears in equation (20) is discussed in the section entitled “Computational Procedure and Discussion.” 
Damping-in-roll coefficient )3 C . — The damping-in-roll coefficient is obtained by proper integration of the pressure 

distribution given by equations (5). The following result is obtained: 
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Rolling-moment-loss coefficient j3 C te . — The rolling-moment-loss coefficient 8 C\ 6 is obtained by suitable integration of 
the pressure distribution given by equations (6) by using 
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from equation (19), since d—y The following result is obtained: 
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Twist parameter 6 r /5 . — The twist parameter 9 r /S can be evaluated from equations (18c) and (18d) for 1—- ~ less than 

or greater than i\ a , respectively. 1& order to perform the evaluation, the value of C Pa is taken from equations (4) to be the 
following: 

c ’-~l7h m 

The use of this linearized-theory result is conservative. The value of k is found from equation (14b) to be the follow ing : 
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COMPUTATIONAL PROCEDURE AND DISCUSSION 

In order to simplify the solution of equations (2) and (3) 
a computational form has been prepared. This computa- 
tional form, which is presented as table I, is practically 
self-explanatory and can be used without reference to the 
analysis. Figures 4 to 15 contain the necessary information 
to be entered in the table. Most of these figures present the 
results of numerical calculations made from equations (20) , 
(21), (23), and (26) ; the variables £ and f appearing in these 
equations have been replaced by the appropriate combina- 
tions of fL4, X, .and c a /c so that the figures contain these 
more familiar quantities. The curves for the rect ang ular 
wing (A= 1.0) were taken from reference 1. Although in 
principle these curves could have been obtained from the 
equations of the present report, the work involved in evalu- 
ating the many indeterminate forms arising when X=1.0 
made such a procedure impractical. 

The factor F which appears in equation (20) is introduced 
to allow for the effect of the t, railing-edge angle. This factor 
is the same one shown in figure 7 of reference 1 and was 
derived by using Busemann’s second-order approximation 
for isentropic flow. For the purposes of the present report, 
the Mach number and trailing-edge angle must be taken 
normal to the aileron hinge line. The factor F is plotted in 
figure 10 against M v , the Mach number component normal 
to the aileron hinge line, with <i> N , the included trailing-edge 


angle measured normal to the aileron hinge line, as the param- 
eter. The relations between M N and M and between 
and i j> are given by the following equations: 


m n — , _ 

Vi+*i 2 



M 

wm 

<t> N = 2 tan -1 1 


—2 tan -1 

1+ ( X 2 e)(li+x) 


(27) 




tan • 


(28) 


For ailerons other than those having flat sides, the trailing- 
edge angle should be taken as the angle between straight lines 
drawn from the trailing edge to the airfoil sin-face at the 
hinge line. For f3C ts this procedure can be shown to give the 
correct value of F independent of the airfoil profile; however, 
this method of determining F for curved-profile contx-ol 
surfaces is not valid for trailing-edge-angle corrections to 
hinge moments. 
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Figure 8.— Continued. 
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• (e) X=0.S. 

(£) X=1.0. (From reference 1.) 
Figtoe 8. — Concluded. 
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Figure 12.— Variation of rolling-moment-loss coefficient with 0A. (From equation (23).) 
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(3A 

(d) X=0.6. (e) X=0.8. (f) X=L0. (From reference 1.) 

Figure 15.— Concluded. 
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The aileron deflection S used in /3(7 {j is measured in a plane 
parallel to the free stream. The value of the rolling-moment- 
effectiveness coefficient referred to an aileron deflection 
measured perpendicular to the hing e line can be obtained, 
if desired, by dividing /3(7 {j by the quantity -y 1-f &i 2 . 

The minimum value of BA for which each curve is plotted 
corresponds in most cases to the limit of applicability of the 
equation for the curve. Extrapolation of the curves is not 
recommended. The maximum value of b a jb for which results 
are given is O.S; this value is expected to lie close to the 
maximum that will be encountered on actual airplanes with 
fuselages. 

The setup of table I can of course be altered. For example, 
a solution may sometimes be required for the torsional stiff- 
ness necessary to maintain a given percentage of the rate of 
roll that would be obtained with a perfectly rigid wing. In 
such a case, the solution can be obtained by recasting the 
last several columns, as shown in the table. If experimental 
data are available for any of the derivatives, these data may 
be used in place of the theoretical values of the present report. 

Although an explicit solution for the reversal speed is 
impossible to obtain because of the complexity of the equa- 
tions, the reversal speed is easily found by calculating the 
rate of roll through the speed range. The variation of the 
rate of roll with speed, however, is usually of more interest 
than the reversal speed alone; therefore, no extra calculations 
are involved in calculating the reversal speed. 

The twist parameter d r /S is directly proportional to q and 
inversely proportional to m» r . (See equations (26).) Be- 
cause of this relationship and because the other derivatives 
appearing in equations (2) and (3) are independent of q and 


m« r , the values of 


P JI2V and _ £ b!2V_ 


found from table I for 


5 (pb(2V ) B ' 
specified values of altitude and stiffness are easily converted 
to values for other altitudes and degrees of stiffness. 

Although the analysis was made by considering two wing 
panels lying in the same plane, experimental evidence exists 
which indicates that the results may with some justification 
be applied without change to configurations consisting of 
three or four wing panels, each containing an aileron, if the 
span and aspect ratio are taken as those of the wing formed 
by the reflection of a single panel. If the configuration con- 
sists of four panels, two of which contain ailerons, then the 
value of fiCi p must be doubled. 

The results of some calculations made for a hypothetical 
configuration are given in figure 16. 

CONCLUDING REMARKS 

An analysis has been presented of the effect of torsional 
flexibility on the rolling characteristics at supersonic speeds 
of tapered unswept wings with partial-span constant-percent- 
chord ailerons extending inboard from the wing tip. 



M 

Figure 16.— Results of computations for a hypothetical configuration at sea level. .4=4; 
X=0.4; ^p=0.4; ^=0-2; 6=30 feet; au, =2,000,000 foot-pounds per radian; 

Because of the complexity of the equations resulting from 
the analysis, numerical calculations were made from the 
equations and the results are presented in a series of figures. 
A computational form was prepared for use with the figures 
so that calculations could be made without reference to the 
analysis. 


Langley Aeronautical Laboratory - , 
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